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Abstract 

Vehicular networks allow vehicles to share information and are expected to be an integral part 
in future intelligent transportation system (ITS). In order to guide and validate the design process, 
analytical expressions of key performance metrics such as packet reception probabilities and throughput 
are necessary, in particular for accident-prone scenarios such as intersections. In this paper, we analyze 
the impact of interference in an intersection scenario with two perpendicular roads using tools from 
stochastic geometry. We present a general procedure to analytically determine the packet reception 
probability and throughput of a selected link, taking into account the geographical clustering of vehicles 
close to the intersection. We consider both Aloha and CSMA MAC protocols, and show how the 
procedure can be used to model different propagation environments of practical relevance. We show 
how different path loss functions and fading distributions can be incorporated in the analysis to model 
propagation conditions typical to both rural and urban intersections. Our results indicate that the 
procedure is general and flexible to deal with a variety of scenarios. Thus, it can serve as a useful design 
tool for communication system engineers, complementing simulations and experiments, to obtain quick 
insights into the network performance. 
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I. Introduction 

Vehicular networks have gained considerable attention in the past years and are regarded as 
one of the key components in future intelligent transportation systems (ITS) [l2l-[|7l- By the use 
of wireless communication they allow vehicles to continuously share information with each other 
and their surrounding (e.g., roadside infrastructure), in order to perceive potentially dangerous 
situations in an extended space and time horizon [l3]|. This enables a new set of applications 
that are expected to enhance both traffic safety and efficiency. These applications include lane 
change assistance, cooperative collision avoidance, emergency vehicle warning, traffic condition 
warning, tolling, hazardous location warning, speed management. 

The IEEE 802.1 Ip standard has been defined to meet the communication demand of these 
applications, and 5G cellular networks standards are being developed to support device-to- 
device (D2D) communication ll8ll- l[T3ll . However, different ITS applications clearly have different 
requirements on the communication links, with the most stringent demands imposed by safety- 
related applications, with extremely low latencies (below 50 ms in pre-crash situations), high 
delivery ratios (for full situational awareness), and relatively long communication ranges (to 
increase the time to react in critical situations) IIT4l - IIT6l . These requirements, in combination 
with a possible high density of vehicles, makes the design of vehicular communication systems 
challenging. This is further exacerbated by high mobility and passing vehicles, which leads to 
rapidly changing signal propagation conditions (including both severe multipath and shadowing) 
and constant topology changes. 

In order to guide and validate the communication system design, extensive simulations and 
measurements are often used m, iini, which are both time consuming and scenario-specific. 
In order to obtain insight in scalability and performance, analytical expressions of key perfor¬ 
mance metrics are necessary, in particular for high velocity scenarios (in particular highways) 
and accident-prone scenarios (e.g., intersections). Stochastic geometry is a tool to obtain such 
expressions, and has been widely used in the design and analysis of wireless networks ifT^ . In 
2-D planar networks, the analysis is well developed and a multitude of approaches to consider 
both geographical and medium access control (MAC) induced clustering [|T9l . Il20l as well as 
different types of fading ll2TI - ll2^ exist. However, in vehicular networks, where the location 
of the nodes are restricted by the roads, previous work that includes the spatial statistics of 
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vehicles typically considers one-dimensional roads Il24l - lf2^ . For these vehicular scenarios, 
geographical clustering has been addressed in ll24ll . while effects due to the 802.lip carrier 
sense multiple access (CSMA) MAC protocol was studied in [[25ll . [|2^ . Il29l . Hence, these 
works enable communication system design for one-dimensional highway scenarios, but do not 
capture well the salient effects of intersections. Intersections were considered explicitly in ifTl , 
[[30l, which found that it is important to properly model the interference from different roads and 
account for the distance of receivers to the intersection, i.e., to take into account the clustering 
of cars around the intersection and the non-stationarity of the spatial distribution. 

In this paper, we present a general procedure for the evaluation of packet reception probability 
and throughput in intersection scenarios, and provide a model repository that can be used to adapt 
to a variety of different environments of importance in the vehicular context. This includes both 
rural and urban scenarios, different propagation conditions, and different MAC protocols. Latency 
and mobility are not treated in this paper. 

The remainder of the paper is organized as follows. Section |II] introduces the system model. 
In Section Hill we discuss typical characteristics of the vehicular channel and show how the 
model can be tailored to different environments. In Section we present a general procedure 
to calculate the packet reception probability near an intersection, as well as the throughput. 
Section rvl shows how the proposed procedure can be used to calculate these performance metrics 
for a number of cases of practical relevance, and how different assumptions on loss function, 
fading, and MAC protocols affect the analytical tractability of the packet reception probability. 
Finally, Section |Vl] summarizes and concludes the paper. 

II. System Model 

We consider an intersection scenario with two perpendicular roads, as shown in Figured] For 
simplicity, we assume that the width of the two roads indicated by H and V can be neglected, and 
that the roads each carry a stream of vehicles, modeled as one-dimensional homogeneous Poisson 
point processes (PPPs). The intensity of vehicles on both roads is denoted by Ah and Ay, and the 
point processes describing the location of the vehicles on the two roads are represented by ~ 
PPP(Ah) and $v ~ PPP(Av). The positions of individual vehicles (also referred to as nodes) 
on the two roads H and V are denoted by x, = [x*, 0]"*^ and x* = [0, yj^ , respectively, assuming 
the roads are aligned with the horizontal and vertical axes. As both vehicle-to-vehicle (V2V) 
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and infrastructure-to-vehicle (I2V) communication is of interest, we eonsider a transmitter (Tx) 
with arbitrary loeation Xtx = The Tx broadeasts with a fixed transmission power 

P. Without loss of generality, we eonsider a reeeiver (Rx) on the H-road at loeation Xrx = 
[xrxjO]^, i.e., at a distanee d = |a:rx| away from the the interseetion ^ The signal propagation 
eomprises power fading S and path loss Z(xtx,Xi.x). At the Rx, the signal is further affeeted by 
white Gaussian noise with noise power N and interferenee from other eoneurrently transmitting 
vehieles on the H- and V-road. The amount of interferenee experieneed by the Rx depends on 
the ehoiee of MAC protoeol. For a given MAC seheme, the position of interfering vehieles at a 
given time ean be represented by the thinned point proeesses and $y We ean express 

the signal-to-interferenee-plus-noise ratio (SINR) as 


SINK 


P S'o/(x^x5 Xj-x) 

^ *S'x^(x 7 ^rx ) + Zlxe-I>“AC PS'x/(x : ^rx ) + N^ 

H V 


( 1 ) 


where Sq denotes the fading on the useful link and S'x denotes the fading on an interfering link 
for an interferer at loeation x. A paeket is eonsidered to be sueeessfully reeeived if the SINR 
exeeeds a threshold (3. 

Our aim is to analytieally eharaeterize (i) the probability that the Rx sueeessfully reeeives 
a paeket sent by the Tx; (ii) the throughput of the link between Tx and Rx. This problem is 
ehallenging due to the speeifie propagation eonditions and interferenee levels experieneed in 
these interseetion seenarios. In the next seetion, we will deseribe these in more detail. 


Remark 1. While the seenario eonsidered here is simple, it ean easily be extended to oases 
where the width of the roads ean not be ignored (without introdueing signifioant modeling 
errors) by splitting the road into several lanes, eaeh modeled as a new road/lane, as disoussed 
in Seetion IV-AI Furthermore, general multi road/lane extensions ean be used to explioitly model 
interferenee from other roads in the surrounding. 


'Note that in the case the Tx is on one of the roads, it can belong to either $h or 4>v (but does not necessarily have to) as 
the results still hold due to Slivnyak’s Theorem (TH Theorem A.5] 

^Note that due to the symmetry of the scenario this also captures the case when Rx is on the V-road 
^For a general MAC scheme, the thinning process is inhomogeneous. 






5 



(a) intersection scenario (b) abstraction 


Figure 1. Illustration of considered scenario: (a) A two-way intersection scenario in which each road carries a stream a vehicles, 
(b) the abstraction used in modeling. The Tx (indicated by the blue car) can be at any location, while the target Rx (green car) 
is located on road H. Other vehicles on the roads H and V, of which some transmit concurrently and cause interference, are 
shown as grey cars. 


III. Models in Vehicular Communication 

Vehicular communication systems must be able to function in a large variety of conditions, 
including in urban canyons and in rural settings. In this section, we will discuss characteristics 
for vehicular channels that are important from an SINK point of view, and detail different models 
regarding path loss, fading, and MAC protocol. 

A. Power decay and blockage 

Extensive measurement campaigns IfT^ . OTTl - OAl have been performed to characterize the 
vehicular channel in a variety of propagation environments such as rural, highway, suburban, 
and urban scenarios. As it is important to understand how the power decays with distance, much 
efforts have been put into finding large-scale path loss models, which characterize the slope of 
distance-dependent power loss in decibels (dB). We will distinguish between line-of-sight (LOS) 
and non-line-of-sight (NLOS) propagation, depending on whether or not the direct LOS signal 
between a Rx and a Tx is blocked. Lor LOS propagation, conventional path loss models, where 
power decays approximately with the squared Euclidean distance between Rx and Tx are well- 
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accepted IfT^ . For NLOS propagation, e.g., in urban eanyons, measurements indieate inereased 
loss over LOS propagation, with eomplex dependeneies on the absolute position of Tx and Rx, 
widths of the roads, and different loss exponents for own and orthogonal road [1^ . If34l . The 
eomplexity of these models render them intraetable when it eomes to mathematieal analysis, 
so we rely on the simpler and more traetable Manhattan model, whieh was first proposed for 
modeling of similar seenarios in the well-known WINNER II projeet [[35]|. In partieular, to allow 
for mathematieal analysis, we suggest the path loss of the following form. 

• For NLOS propagation, where the direet line-of-sight (LOS) between the Rx and the Tx is 
bloeked by buildings and the signals have to propagate along the urban eanyons formed by 
the orthogonal streets, we use the Manhattan distanee: 

(^txj ^rx) ^ ll^rx ^tx|li ; (2) 

where ||xrx — Xtx||i is the norm, a > 0 is the path loss exponent, and A is a eonstant that 
depends on several faetors sueh as antenna eharaeteristies, earner frequeney, and propagation 
environment. 

• For LOS propagation, where the direet LOS between the Rx and the Tx is unobstrueted, 
we use the Euelidean distanee 

^E(^tX5 ^rx) ^ ll^rx ^tx||2 ; (3) 

where Hx^x — Xtx ||2 is the ^2 norm. Note that the values of a and A might be different in 
dH) and dH). 

B. Random power variations due to fading 

Fading refers to random fluetuations in the reeeived power around the average reeeived 
power, given by the path loss. The fading experieneed on a link depends on the seenario 
and the environment, and is typieally modeled as a random variable If36ll . For example, near 
a rural interseetion, vehieles are likely to eommunieate via LOS links, and exponential fading 
is eonsidered an appropriate model [[8]|, lf34l . On the eontrary, if the interseetion is loeated in 
an urban environment with tall buildings, the fading for NLOS links is modeled using a log¬ 
normal model [l33l|, lf34l . with typieal values on power variations with respeet to the path loss for 
NLOS interseetions are in the range of 3-6 dB. Based on these empirieal results, we will model 
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the fading as log-normal (and approximated by an Erlang random variable for mathematical 
tractability - see Section ITV-A2I) for NLOS links and as exponential for LOS links. 

C. MAC protocols 

The MAC protocol governs when a user can access the channel, and aims to control the 
interference in the network. The two most common MAC protocols for ad-hoc networks are 
Aloha and CSMA. In Aloha, which is the simpler of the two, nodes that have a packet to send, 
access the channel during a time slot with a probability p G [0,1]. In contrast, in CSMA, before 
sending a packet, a node verifies that the channel is free by listening to the channel. Only if 
the channel is free, the node transmits the packet. If the channel is busy, the node is forced 
to wait a random back-off time before it can try again lUTll . The 802.lip standard, which has 
been designed for the first generation vehicular networks, will rely on a CSMA/CA (collision 
avoidance) MAC. We will consider CSMA as well as Aloha, as Aloha is easier to analyze and 
has been argued to exhibit similar performance as CSMA, especially for dense networks ll26l . 

IV. Stochastic Geometry Analysis 

From Section Hill it is apparent that vehicular communication systems will operate under a 
variety of propagation conditions. In this section, we describe a general and unified methodology 
to compute the communication performance for all these conditions, as well as different MAC 
protocols. In particular, we will determine (i) the packet reception probability P(/9, x^x, Xtx), 
i.e., the probability that a receiver located at x^x can successfully decode a transmission from a 
transmitter located at Xtx, in the presence of interferers on the H- and V-road; (ii) the throughput 
T(/3, Xj-x, Xtx), i-C., the expected rate for the link between the Rx and Tx at locations x^x and Xtx, 
accounting for both the packet reception probability and the probability of gaining access to the 
channel. Both P(/3, Xj-x, Xtx) and T(/9, x^x, Xtx) depend on the loss function, fading distribution, 
and the MAC protocol. Note that the loss function and fading distribution relate to the power 
decay and blockage as well as the random signal variations in the specific scenario, while the 
MAC protocol relates to number of interferers and their locations. Several applications of this 
methodology will be discussed in Section |Vl 




A. Packet reception probability 


To derive the paeket reeeption probability for the interseetion scenario, we start by accounting 
for the fading distribution of the useful link. We express 

P(Ax,.,xt.) =Pr(SINR>/3) 

= Pr (^0 > (/h + Iy + N) (4) 

in which N = N/P and Jh = S'x/(x, x^x) while Jy = 5'x/(x, x^x). Condi- 

H V 

tioning on the path loss, we can now write the packet reception probability as 

P(/?,X,x,Xtx) (5) 

= |-^5o + ly + n'^ /?//(xtx, Xj-x) j I 

= JJ ^So + ^2 + //H,rv(^i5 f2)dfidf2, 

where (3 = /3//(xtx, Xrx) and Fs^{so) is the complementary cumulative distribution function (CCDF) 
of the random variable S'o, evaluated in sq- 

The expression dS]) can be interpreted in two ways: (i) as the expectation of Fso{{Ih + ly + 
N)f3/li Xtx, Xrx)) with respect to the interference distribution; and (ii) as the transformation of the 
interference distribution with a kernel function determined by the CCDF of the fading distribution 
of the useful link. In either interpretation, the distributions of the interference and the fading 
play an important role. Note that for all relevant fading distributions of the useful link, ([5]) will 
result in the Laplace transform (LT) of the interference distribution or a function of LTs of the 
interference distribution. It is therefore convenient to express these distributions through their 
(LT) or, equivalently, their moment generating function (MGF). 

1) LT of the interference: From (jS]), we see that the packet reception probability P(/9, Xj-x, Xtx) 
is a function of the interference distribution, which itself depends on the location of the Rx and 
the interferers, as well as their fading distributions and path loss. For a general MAC protocol 
the interference from the H- and V-road are not independent. However, for the MAC protocols 
studied in this paper the interference distribution factorizes as fi^jy{ti,t 2 ) = fi^{ti)fiy{t 2 ). In 
fact, the interference is independently thinned on the H- and V-road in the case of Aloha, while 
for the CSMA scheme we can approximate the joint interference distribution as the product of 
the marginals, where the dependence is captured by a location dependent thinning of the original 


PPPs IfJTlI . This means that the interfering point proeesses and either are, or are 

approximated as PPPs (for more details see Seetion IIV-A3I) . and that we ean foeus on a single 
road R G {H, V}, with interferenee distribution fj^. The Laplaee transform of fj-^ is defined as 


^/r(s) = E[exp(-s/R)], 


( 6 ) 


in whieh 

Ir= ^ S'x((x,Xrx). (7) 

^g^MAC 

Substitution of dV]) into dH) then yields 




— Ed 


JJ Es^ {exp (-S S'x/(X, Xrx))} 


^g(j,MAC 


n >C5^(s^(x,Xrx)) 


(b) 


X^^MAC 
/•+00 


= exp (- / (x (z) , Xtx) (1 - Cs^ (s i(x(z),x,^))) dz j , 


( 8 ) 

(9) 

( 10 ) 


where (a) holds due to the independenee of the fading parameters, Ed> [•] is the expeetation 
operator with respeet to the loeation of the interferers, and £ 5 ^ (•) is the LT of the fading 
distribution of the interfering link; (b) is due to the probability generating funetional (PGFL) for 
a PPP [R^ Definition A.5], in whieh (^(^)) ^tx) represents the intensity of the PPP 

whieh depends on the speeifie MAC protoeol and in some eases on the transmitter’s loeation. 
Note that in (fTOl) . the intensity is defined over 2 ; G M, whieh represents the position along the 
road R G {H, V }, where 

x(z) = 

To determine £f^(s), we must be able to eompute the integral (fTOl) . whieh involves knowledge 

of A^^^ (x (z) , Xtx) and £ 5 ^ (s). 


[zO]T ,R = H 
[Oz]T ,R = V 


( 11 ) 


Remark 2. The Laplaee transform of the interferenee ean also be eomputed using the prineiple 
of stoehastie equivalenee [|231 . where the LT in ease of an arbitrary fading distribution ean be 
found based on the LT in ease of Rayleigh fading, given an appropriate sealing of the system 
parameters. 
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2) LT of fading: For many relevant fading distributions, the LT is known, ineluding for 
exponential. Gamma, Erlang, and random variables. While the log-normal distribution is 
harder to deal with, it ean be approximated by the Erlang distribution Il38]| . whieh eombines 
traetability with expressiveness. When S'x ~ E(fc,6'), i.e., an Erlang distribution with shape 
parameter /c G N and rate parameter 1/9 > 0, then 

>^^ 5 . (s) = (1 + sOf . (12) 


As speeial eases, (i) fc = 1 eorresponds to an exponential distribution with mean 6*; (ii) 6 = 1/k 
eorresponds to Nakagami-m power fading. 

Remark 3. When the fading of the useful link is exponentially distributed, dS]) allows us to inter¬ 
pret P(/9, Xrx, Xtx) as the ET of the interferenee, so that P(/9, Xrx, Xtx)= exp [-N~^/e)ci,0/e)Ci,0/e). 

3) Intensity of the interfering PPPs: The intensity (x ( 2 :) ,xtx) of the interferenee de¬ 
pends on the type of MAC that is utilized. We distinguish between two oases: Aloha with transmit 
probability p G [0,1], and CSMA with oontention region with radius 5 > 0. 

• Aloha: Eor an Aloha MAC, the vehioles on eaeh road will transmit with a probability 
p. This leads to an independent thinning of the PPPs, so that (x (z ), Xtx) = pAr, 
irrespeotive of z or xtx. 

• CSMA: Eor a CSMA MAC, a vehiole will transmit if it has the lowest random timer within 
its sensing range (oontention region). This means that (i) the intensity is in this ease also 
a funetion of Xtx as other nodes in its oontention region are foreed to be silent when it is 
aetive; (ii) the interferenee from the H- and V-road is not independent. The timer proeess and 
the oorresponding dependent thinning result in a Matern hard-oore proeess type II, whieh 
ean be approximated by a PPP with independently thinned node density. The approximation 


of the hard-eore proeess by a PPP is shown to be aeeurate in 
the oontext of heterogeneous oellular networks, for instanoe in 


and has been applied in 
When the transmitter 


at Xtx is aetive the resulting intensity of the PPPs used to approximate the point proeess of 


''The extension to CSMA schemes with discrete back-off timers has been proposed in 125II . which retains concurrent transmitters 
due to the non-zero probability of nodes with the same timer value. 
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interferers can be expressed as 



Pa(x(2:))Ar ||x(2:) 


0 


x(z) -Xtxll < s 


- Xtx 


> 5 


(13) 


In (flBl) . Pa (x (z)) is the access probability of a node. The access probability (which is used 
to thin the original process) is the probability that the given node has the smallest random 
timer in the corresponding contention region (in this case modeled as a 2-dimensional ball 
B 2 {^{z), 6 ) with radius 5 centered at location x (^)), and can for one of the roads be 
expressed as 



(14) 


1 -exp(-A(g2(x(z) ,(5))) 
A{B2{^{z),6)) 


(15) 


where 



(16) 


represents the average number of nodes in the contention region. Note that the average 
number of nodes, and thus the access probability depends on the position 2 ; along the 
road and the intensities Ar and Ar/, which here represent the intensities of the unthinned 
processes on the relevant road R and the other road, respectively. 

B. Throughput 

From a system perspective, the packet reception probability is not sufficient to characterize 
the performance, since a MAC that allows few concurrent transmissions leads to high packet 
reception probabilities but low throughputs. Thus, to be able to compare the impact of different 
MAC protocols, we characterize the throughput for the intersection scenario, i.e., the number 
of bits transmitted per unit time and bandwidth on a specific link. For the general case with a 
receiver and transmitter located at x^x and Xtx, respectively, we express the throughput as 


T (/), Xrx, Xtx) = PA(Xtx)P(/3, Xrx, Xtx) loga (1 + /?) 


(17) 


where Pyi(xtx) is the access probability of a transmitter located at Xtx, i.e., the probability that 
the transmitter obtains access to the channel to transmit a packet. For the Aloha MAC, the access 
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probability is simply p^(xtx) = p, while for the CSMA ease the aeeess probability is given in 
(fTSl) and depends on the void probability in the 2 -dimensional ball used to model the eontention 
region around Xtx- 

C. General Procedure 

Given the analysis in the previous subseetions, the general proeedure for determining the 
paeket reeeption probability P(/3, Xrx, Xtx) and the throughput T(/9, Xrx, Xtx) is thus as follows: 

• Step 1: Determine the fading LT £ 5 ^ (s) for the interfering links, as deseribed in See- 
tion IIV-A2I 

• Step 2: Determine the intensity of the interferenee PPP (x (z) , xtx) for R G {H, V}, 
as deseribed in Seetion IIV-A3[ 

• Step 3: From step 1 and step 2, determine the LT of the interferenee for R G {H, V} 

using (fTOl) . 

• Step 4: Determine the fading LT Csq (s) for the useful link, as deseribed in Seetion irV-A21 

• Step 5: From step 4 and step 3, determine P(/3, Xrx, xtx) using ([5]), either by drawing samples 
from the interferenee, or by eonsidering the CCDF of the fading on the useful link as a 
kernel in a transformation (i.e., evaluating a funetion of LTs of the interferenee distribution). 
Finally, use the obtained paeket reeeption probability P(/9, Xrx, xtx) in eonjunetion with the 
aeeess probability p^(xtx) used in step 2 to determine the throughput T(/9, Xrx, Xtx). 

Whether or not eaeh step is traetable depends on the assumptions we make regarding the loss 
funetion, the fading distribution, and the MAC protoeol, whieh will be further diseussed in 
Seetion |Vl 


V. Case Studies 

In this Seetion we present three ease studies to show how the different models presented in the 
paper ean be used to model both rural and urban interseetion seenarios, and how shadowing, LOS 
bloekage, and different MAC protoeols affeet the performanee of the eommunieation system. In 
Case I, we present the most basie ease whieh eorresponds to the rural setting, while in Case II, 
we show how an urban interseetion ean be modeled. Finally, in Case III, we will study the 
impaet of the different MAC protoeols. In eaeh ease study, we will diseuss the traetability of the 
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resulting expressions, validate modeling assumptions through simulations, and provide numerieal 
performanee examples. 

A. Case I - Rural intersection with Aloha 

In the rural interseetion seenario |[I1, [1^ . vehieles are assumed to eommunieate via LOS 
links. Henee, path loss is deseribed by the Euelidean distanee loss funetion /e(')’ defined in ([3]) 
with path loss exponent a = 2, while power fading is modeled with an exponential distribution 
(i.e., S ~ £'[1,1]), for both useful and interfering links. Furthermore, we eonsider an Aloha 
MAC with transmit probability p. 

1) Packet reception probability : Using the proeedure from Seetion ITV-Cl the paeket reeeption 
probability for the rural interseetion seenario is given in Proposition |4] (see also [|T1|, [[30ll l. 

Proposition 4. Given a slotted Aloha MAC with transmit probability p, exponential fading (i.e, 
S ~ £(1, 1)) for each link, Euclidean loss function l-^f) with path loss exponent a = 2, and a 
scenario as outlined in Section m the packet reception probability can be expressed as 



(18) 


Proof: See Appendix 


We note that the paeket reeeption probability eomprises three faetors: the first faetor eorre- 
sponds to the paeket reeeption probability in the absenee of interferers; the seeond faetor eaptures 
the reduetion of the paeket reeeption probability due to interferers on the H-road; the third faetor 
eaptures the additional reduetion of paeket reeeption probability due to interferers on the V-road. 

Remark 5. Proposition |4] ean be extended in a number of ways: 

• As was noted in [[T]| additional roads/lanes with arbitrary orientations ean be aeeounted for, 
eaeh road eontributing with an additional faetor to the paeket reeeption probability. This 
approaeh ean for example be used to take into aeeount interferenee from surrounding roads. 
Furthermore, it ean be used to handle eases where the width of the roads ean no longer be 
ignored, by splitting the road into several lanes. 
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• Extensions to scenarios with non-homogeneous PPPs are also possible, in order to model, 
e.g., clustering of vehicles due to traffic congestions. In general this requires numeric 
integration to evaluate the LTs of the interference distribution, but for special cases such as 
piecewise linear intensity functions, closed-form expressions can be found. 

2) Numerical results: Throughout the remainder of the paper, we consider an intersection 
where the intensity of vehicles on the two roads are Ah = Ay = 0.01 (i.e., with an average 
inter-vehicle distance of 100 m). Furthermore, we assume a noise power N of —99 dBm, 
an SINR threshold oi (3 = 8 dB [fTTl . and that A = ?> ■ 10“®, approximately matching the 
conditions in ll^ . We set the transmit power to P = 100 mW, corresponding to 20 dBm. For the 
purpose of visualization, we show the outage probability Pout(/9, Xj-x, Xtx) = 1 — P(/3, Xrx, Xtx) 
instead of the packet reception probability. Figure shows the analytical outage probability 
for the rural intersection scenario as a function of distance between transmitter and receiver 
||xrx — Xtx ||2 for different distances to the intersection d G {0 m, 100 m, 500 m} and different 
transmit probabilities p G {0, 0.005, 0.1}. We observe that the outage probability increases with 
the distance between the receiver and the transmitter, and that interference has a negative impact 
on the performance as the outage probability is higher for increased transmit probabilities. In 
the absence of interferers (p = 0 ) the system achieves an outage probability of 10 % when 
the receiver and transmitter are spaced approximately 600 m apart. When p is increased to 
0.005 the communication range is drastically reduced to about 130 m, due to the interference. 
Furthermore, the figure reflects the location dependence of the outage probability with respect 
to the intersection, and we can see that the outage probability increases when the receiver is 
closer to the intersection. For the purpose of validation, we have added Monte Carlo simulation 
with 10,000 realizations of the PPPs and fading parameters, perfectly matching the analytical 
expressions. 

B. Case II - Urban Intersection with Aloha 

This case study, model an urban intersection scenario, with the Tx on the V-road and the Rx 
on the H-road. Signals arriving to the Rx from the V-road are assumed to be in NFOS, modeled 
through Manhattan path loss and Erlang fading (which serves as an approximation of log-normal 
fading). Signals arriving to the Rx from the own H-road are in EOS, modeled through Euclidean 
path loss and exponential fading. 
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Figure 2. Comparison of analytical and simulated outage probability Pout(/3, x^x, Xtx) versus distance between transmitter 
and receiver ||xrx — Xtx ||2 for different distances to the intersection d as well as different transmit probabilities p. The distances 
d are 0 m (diamonds), 100 m (circles) and 500 m (squares). 


1) Packet reception probability: The packet reception probability for the urban intersection 
scenario is given in Proposition [6l 

Proposition 6. Given a slotted Aloha MAC with transmit probability p, Erlang fading (i.e., 
S ~ E{ko,9o)) and Manhattan loss function /m(‘) for the useful link, Erlang fading (i.e., S ~ 
E{kv,0v)) and Manhattan loss function for the interfering links from the V-road, exponential 
fading (i.e, S ~ E(l, l)j and Euclidean loss function /e(') for the interfering links from the Id- 
road, and a scenario as outlined in Section m the packet reception probability can be expressed 


as 


/cq 1 ^ / ' \ /"i 


i=0 j=0 


(19) 


where 
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a 
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^v, —l-Q', IH I" ~ 


/r which K = 2pAHx4^/“7r/Q! CSC (vr/a) and C = /9 ||xrx — Xtx||i / (x46'o) . 


Proof: See Appendix E ■ 

We observe that the analytical expressions become more involved when changing the loss 
function as well as the fading distribution for the links to the V-road, but in contrast to the 
rural intersection scenario it is possible to obtain closed form expressions for a general a (this 
is because Manhattan path loss for the interferers from the V-road is easier to handle than 
Euclidean path loss). Furthermore, it should be noted that if the Tx is assumed to be on the 
H-road, the expressions become more compact. Moreover, similarly as for the model presented 
in |[34l . Proposition only gives realistic results when the Rx and the Tx are at least a few 
meters away from the intersection. This is because when the Rx is at the intersection, all links 
become LOS, while when the Tx is at the intersection, the useful link becomes LOS. In either 
case, the corresponding links should be modeled with exponential fading, rather than Erlang 
fading. 

2) Numerical results: In this section we intend to validate the accuracy of the Erlang approxi¬ 
mation. We consider the same parameters for the LOS propagation as in Section IV-A[ We set the 
Aloha transmit probability to p G {0.002, 0.02}. For all NLOS links, we use the same value of 
A as in the LOS links, set a = 2 , and consider the fading to be log-normal with 3.2 dB standard 
deviation, as in one of the intersections studied in lf33l . Maximum likelihood fitting of the Erlang 
distribution to the log-normal distribution yielded ko = ky = k = 2 and Oq = 9y = 9 Ri 0.66. In 
Figure m a comparison between the outage probability obtained by evaluating Proposition [funder 
the Erlang approximation, and Monte Carlo simulations for the same scenario but with 3.2 dB 
log-normal shadowing is shown. First, we observe that the analytical results based on the Erlang 
approximation agrees well with the simulations, i.e., in terms of outage probability the Erlang 
fading provides a good approximation to the actual log-normal fading. Furthermore, we see that 
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Figure 3. Comparison between analytical and simulated outage probability Pout(/?, Xrx, Xtx) versus distance between receiver 
and intersection d for different Tx positions and transmit probabilities. Lines show analytical results based on Proposition 
while markers show simulations for the same scenario but with 3.2 dB log-normal fading. The different transmit probabilities p, 
which the outage probability is plotted for are 0.002 (solid) and 0.02 (dashed) lines. 


as expected, lower transmit probability, as well as placing the Tx closer to the intersection results 
in lower outage probabilities. Even though the results shows very good agreement between the 
analytical results and the simulations, it should be mentioned that the approximation become less 
accurate when the standard deviation increases. In particular, this is noticeable when the standard 
deviation of the log-normal distribution exceeds 3.8 dB, as the Erlang distribution obtained from 
the fitting then reverts to the exponential distribution. Eurthermore, when decreasing the standard 
deviation of the log-normal distribution, the estimated value of k rapidly increases, so that the 
model becomes less tractable. 

C. Case III - Aloha vs CSMA 

In this final case study, we will focus on the MAC protocol and how it affects performance 
and tractability. To do this, we start from the rural intersection scenario, but replace the Aloha 
MAC with a CSMA MAC. As the MAC affects not only the packet reception probability but 
also the access probability, we will also consider throughput in this case study. 
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1) Packet reception probability: The packet reception probability for the CSMA case is given 
in Proposition |71 


Proposition 7. Given a CSMA MAC with contention radius 6, exponential fading (i.e, S ~ 
E{1,1)) for each link, Euclidean loss function /e(‘) with path loss exponent a = 2, and a 
scenario as outlined in Section m the success probability can be expressed as 

P(/3,Xrx,Xtx) = e~~Cj^0)Ci,,{/3), (22) 


where (3 = /9/ZE(xtx, Xrx), and 


£/jj(s) = exp 


'•+0O ,xtx) 

-oo 1 + lairx - xf /As 



(23) 


Ciy{s) = exp - 


r.+oo (^[0,2/]'^ ,Xtx) 

-oo 1 + ||[Xrx, 


dy 


(24) 


where 0]""^, Xtxj and Ay^*" ^[0, y]"'", Xtxj are given in and ^IbTl) . respectively. 


Proof: See Appendix O ■ 

As can be seen from Proposition |71 the expressions we obtain still involve an integral that 
can be solved numerically easily and efficiently. 

2) Numerical results: In order to evaluate the accuracy of the approximation introduced 
in Section IIV-A31 we start by comparing the analytically calculated outage probability to a 
simulation with 50000 realizations of the fading parameters and the hard-core process induced 
by the dependent thinning resulting from the CSMA scheme. This comparison can be seen in 
Figure S which shows the analytical and simulated outage probability as a function of the dis¬ 
tance between the receiver and the intersection for two different transmitter locations (xtx = [0,0] 
and Xtx = [0,150]), as well as two different CSMA contention radiuses 5 G {500 m, 10000 m}. 
We observe good correspondence between simulation and analytical results, and an increase in 
outage probability with increased distance to the intersection. We also note that when Xtx = [0,0], 
it is possible to compare Figure |4| with Figure [2l We note that for 5 = 10000 m, for a distance 
of 100 m between Rx and intersection, CSMA has an outage probability of 0.003, while Aloha 
is over 25 times worse, with an outage probability of 0.08. 
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distance between receiver and intersection, d [m] 

Figure 4. Comparison of analytical (blue lines) and simulated (red markers) outage probability Pout(/^, Xrx, Xtx) versus 
distance between receiver and intersection d, for different transmitter locations Xtx as well as CSMA contention radiuses 5. The 
receiver is located on the Fl-road, while the transmitter location is fixed to either xtx = [0, 0] (solid lines) or xtx = [0, 150] 
(dashed lines). The different CSMA contention radiuses are S = 500 m and S = 10000 m, which in the region where the access 
probability is constant, i.e., far away from the intersection, corresponds to pA = 0.1 and pA = 0.005, respectively. 


To further study the performanee gains aehieved by using CSMA eompared to Aloha, we 
now look at both outage probability and throughput for a speeifie reeeiver and transmitter 
eonfiguration. The eonfiguration that we eonsider is Xrx = [0 0]"’" and xtx = [Rm mm O]"*". Note that 
for the Aloha ease this plaeement results in the worst possible throughput for a fixed x^x)- 

Figure [5] and Figure [ 6 l show the outage probability as well as throughput as a funetion of the 
aeeess probability pAi^tx), for two different values on i?comm ^ {100 m, 200 m}. 

For Aloha (Figure [5]), we see that with an inerease in pA^Xtx), outage probability inereases 
due to the presenee of more interferers. The throughput first inereases (due to more aetive 
transmitters) and then deereases (due to overwhelming amounts of interferenee), leading to an 
optimal value of p^(xtx). However, in order to guarantee a eertain quality of serviee, one must 
also eonsider a guarantee on the outage probability. For instanee, if we want to guarantee an 
outage probability of less than 10 % on the link when i?comm = 100 m, the optimal value of 
PA(xtx) ~ 0.006, leading to a throughput of around 0.0055 bits per unit time and bandwidth. 

For CSMA (Figure[ 6 l), a low aeeess probability (i.e., large eontention region) reduees the outage 
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Figure 5. Outage probability Pout (/3, Xrx, Xtx) and throughput T {p, Xrx, Xtx) for the Aloha case as a function of the transmitter 
access probability pA(xtx). The receiver is located at Xrx = [0, 0], and solid lines correspond to Rcomm = 100 m, while dashed 
lines correspond to Rcomm ~ 200 m. The red circles indicate the maximum throughput that is possible to achieve while 
guaranteeing that the outage probability is kept below the target value of 10 %. 


probability. Similar to Aloha, the throughput first increases with increased access probability and 
then decreases. To achieve an outage probability below 10 % when i?comm = 100 m, the optimal 
value of pA(xtx) ~ 0.023 (corresponding to a contention radius 6 of about 1100 m), results in 
a throughput of about 0.059 bits per unit time and bandwidth. Hence, in this scenario, using 
CSMA instead of Aloha leads to more than a tenfold increase in the throughput for the same 
communication range. 


VI. Conclusions 

We have provided an overview of the dominant propagation properties of vehicular communi¬ 
cation systems near intersections, for both rural and urban scenarios. Based on these properties, 
we proposed a general procedure to analytically determine packet reception probabilities of indi¬ 
vidual transmissions as well throughput, mainly applicable to 802.1 Ip communication. In contrast 
to traditional cellular networks, the one-dimensional road geometry leads to non-homogeneous 
packet reception probabilities and throughputs. 

We have applied and validated this procedure to three case studies, relevant for vehicular 
applications. The results indicate that the procedure is sufficiently general and flexible to deal 
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Figure 6. Outage probability Pout(/?, Xrx, Xtx) and throughput Xrx, Xtx) for the CSMA case as a function of the 

transmitter access probability pA(xtx). The receiver is located at Xrx = [—100, 0], and solid lines correspond to -Rcomm = 100 m, 
while dashed lines correspond to -Rcomm = 200 m. The red circles indicate the maximum throughput that is possible to achieve 
while guaranteeing that the outage probability is kept below the target outage probability of 10 %. 


with a variety of scenarios, that its performance results match well with simulations, and that 
it is able to capture the clustering of vehicles close to the intersection. The procedure can thus 
serve as a useful design tool for communication system engineers, complementing simulations 
and experiments. 

Possible avenues for future research includes validation of the model agains actual measure¬ 
ments, adoption of advanced MAC schemes as well as 5G D2D features. 

Appendix A 

Proof of Proposition 4 

In order to determine the packet reception probability when S ~ i?(l, 1), we follow the 
general procedure from Section llV-Ci 

Step 1: The fading LT for the interfering links can be expressed as Cs^ (s) = 1/(1 + s). 

Step 2: According to Section ITV-A3I the intensity of the two interfering PPPs and 
are ([xO]"^, Xtx) = rAr and A^^^^ ([0, Xt^) = pAy, respectively. 

Step 3: The LT of the interference for the two roads are derived as follows. Since the fading 
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of the interfering links is exponentially distributed, (fTOl) simplifies to 


Ci^{s) = exp (^- J 




-dz 


(25) 


1 + 1 /(s 6 '/(x( 2 :) ,x„)) 

Using (1251) for the horizontal road with Euelidean path loss, and bearing in mind that x(z) = 
[ 2 : 0 ]"'" , we ean write 


( n+oo 

- / 

J —OD 


= exp ( -pXn 


— OD 
+00 


1 + llXrx 
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XU, 


(b) 

= exp 


- 2 pAh 


1 + l^^rx — X\°' /As 
+ 00 


rjAs 

dx 


dz 


1 + 


du 


(26) 

(27) 

(28) 
(29) 

where (a) uses the faet that the intensity is pAn on the H-road, and (6) involves a ehange of 
variable u = |xrx — x\ /{AsY^^. For the partieular ease of a = 2 the FT of the interferenee 


= exp ( —2 pAh (As)^^“ vr/a CSC (vr/a) 


further simplifies to 


For the V-road, using (1251) . we ean in a similar way as for the H-road write 


F/h(s) = exp (-pAhvt-v/As 


(30) 
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Ci^{s) = exp 

Now using that the distanee 


'‘+00 




-pAv 
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1 + 
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-y] 2 Ms 


-dp 


(31) 
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[xrx, —y]^ = y/ x"/^ + p2 = \/d? + we ean introduee Vy = 

d? + p^, with dvy/dy = y/vy. Noting that a PPP remains PPP under a non-linear transformation 
aeeording to the mapping theorem ifTSl Theorem A.l], we have 


F/v(s) = exp^2pAv 
i(^pAv (Als)^/" 


r*-\-00 


= exp 
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(33) 


V ' ■ ' ' Jcoo V^-^0 (1 + 

where we have earried out the following ehange of variable u = {vy/ (y4s)^^“j , and further in- 

trodueed cuo = {d/ (y4s)^^“ j . For a = 2, the integral ean be eomputed as /^°°(v'a; — cuo (1 + ce))' 
vr/y/r+Tu, whieh yields 


Mce 
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Note that for —)■ 0, (1^ reverts to (1^ . while for d —)■ +cx), (1^ tends to one. 

Step 4: The fading on the useful link is eharaeterized by its LT Cs^ (s) = 1/(1 + s) and 
CCDF Fsoi^) = exp (—s). 

Step 5: Using the LT of the interferenee from Step 3, and the CCDF of the fading from 
Step 4 as a kernel, we ean now determine P(/3, Xrx, Xtx) through ([5]). First using the CCDF, and 
evaluating it in the desired point, we ean write 

^So + ^2 + = exp (^ti + ^2 + (35) 

As the interferenee from the H- and V-road is independent (i.e., and are independent) 
we ean now use (1351) to express the transform in ([5]) as 

P(/3, Xrx, Xtx) 

/*H”CXD /*-1-C50 

= exp (-iV/?) J exp (-ii/s) dti J fiyih) exp (-* 2 / 3 ) di 2 

= {J}^ 

Using the results from step 3, and the variable ehanges /3 = /)//E(xtx, Xrx) and N 
finally allow us to express the paeket reeeption probability as ([18]). 

Appendix B 

Proof of Proposition 6 

We use the proeedure from Seetion IIV-CI 
Step 1: The fading LTs for the interfering links from the H-road and the V-road ean be 
expressed as £ 5 ^ (s) = 1/(1 + s) and £ 5 ^( 5 ) = 1/(1 + sOvY^ , respectively. 

Step 2: According to Section HV-A3I the intensity of the two PPPs and <hyare pXu 

and pAv, respectively. 

Step 3: The LT of the interference for the two roads are derived in the following way. For the 
H-road, with interferers x ’ *^he fading LT as well as the loss function are the same as 

in the rural intersection case. Thus we can according to (l29l) expresse the LT of the interference 
for a general a as 

£/h(s) = exp 2 pAh (^s)^'^“ 7 r/a CSC (vr/a) j . (38) 

For the V-road we now have fading LT £s„(s) = 1/(1 -f s9yY^ , intensity pAy, and Manhattan 
loss function. Hence, using ([TO we can write 


(36) 

(37) 

= N/P, 
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(39) 

(40) 

(41) 

(42) 

(43) 


where (a) uses the Binomial Theorem and the variable change s9yA —)■ b, ( 6 ) uses that for 
points X G the distance ||xrx —x||j^ = |a;rx| + |p| = d + \y\, and (c) uses the variable 

change d + |p| —)■ u. For q>0,ky>q+l, b>0 and d > 0 the integral can be evaluated in 
closed form, and for a general a we can express the LT of the interference as 


d^ivY) = 
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where 2-^1 is the regularized hypergeometric function. Note that for a = 2 and ky = 6y = 1 
(i.e., exponential fading) this simplifies to 


Ci^{s) = exp pAv\/^ ^TT — 2 arctan^—^^ , 


(45) 


and when d —)■ 0 we get 


/)/v(<s) = exp {^—yXyn^fA^ , 


(46) 
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i.e., it reverts to the same form as Ci^{s) in ( 1 ^ . 

Step 4: The fading on the useful link is eharaeterized by its LT Cs^ (s) = 1/(1 + s6o)^° and 
CCDF 


fco —1 


= e-/"" y; - 


i=0 


ilOl 


(47) 


Step 5: In the same manner as in Appendix we now use the LTs of the interferenee from 
Step 3, and the CCDF of the fading from Step 4 to determine P(/), x^x, Xtx) through ([5]). First 
using the CCDF, and evaluating it in the desired point, we ean write 


F, 


*50 


= e 


+ F + Nj /) 

fco —1 


1 


E Titf ('’) {F+h+N 

i=0 0 


(a) _ 

= e 


ko-l 

C(4l+t2+Af) EM h+t2+N 


i=0 


g-C^ 


fc 0 ~l 7 


eh;) i\ ( i) 2 5 


(48) 

(49) 

(50) 

(51) 


i=0 j=0 

where (a) involves the variable ehange ( = 13/9q and ( 6 ) uses the Binomial Theorem. Due to 
the independenee of the interferenee we ean now use (I5TI) to express the transform in (|5]) as 

ko — l i 
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are obtained using the Laplaee transform property t"/ (t) <—>■ (—1)" [f (f)] (C)- Note that 
(l52l) and (1551) use the variable ehange N = N/P. Now using the results from step 4 express the 

^th 

derivative of the LT of the interferenee from the H-road as 


^c, to = V V 

d{C)^ « m\{-m + l)\ 


(59) 


where (■)^ is the Poehhammer symbol and k, = 2pX}i tt/o esc (vr/a). For the V-road, there 

is no general eompaet expression for the derivative of but an explieit expression ean 

in prineiple be ealeulated for any n, ky and 9y. 


Appendix C 

Proof of Proposition 7 
We use the proeedure from Seetion IIV-CI 

Step 1: The fading LT for the interfering links ean be expressed as Cs^ (s) = 1/(1 + s). 

Step 2: Aeeording to Seetion HV-A31 the intensity of the two PPPs <hy and are for 

this ease also a funetion of the transmitter loeation x^x- Using (fT3]) we ean express the intensity 
for the H-road as 

l-exp(-2^AH) -y c T? 

25 ^ '''1 

Ar" ([X. 01-, x„) = X e K. (60) 

0 else 
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in whieh TZi = {a;| |a;| > 6 and \/(x — XtP)"^ + Vtx > <5} and 7^.2 = {a;| |a;| < 6 and \/{x — Xt^y + i/tx > 

(5}. Similarly for the V-road, 
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6 }. 

Step 3: Using (1251) the LT of the interferenee for the H- and V- road ean be expresed as 
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and 


Ci^{s) = exp 


.+00 (^[0,7/]^,Xtx) 
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1 + 


[Xrx, -yf /As 


-dy 


(63) 
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Step 4: The fading fading on the useful link is eharaeterized by its LT Cs^ (s) = 1/(1 + s) 
and CCDF Fs^is) = exp (—s). 

Step 5: By applying a loeation dependent thinning, we approximate the interferenee from 
the H- and V-road as independent. Henee, as the fading on the useful link is exponential (i.e., 
So ~ E{1, 1)), we ean in the same way as in Appendix lAl express the paeket reeeption probability 
as P(/3, Xrx, Xtx) = e' -fKc,, (/3) C,, ( /3 j. Using the results from Step 3, as well as the variable 
ehange N = N/P, we can for the particular value of a = 2 finally obtain (l22l) . Note that for 
a general transmitter location Xt^, we are not able to evaluate the integrals in (1^ and (1^ in 
closed form, but have to resort to numerical evaluation. 
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